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^ Abstract 
> 

We propose to use the recoilless (d,^He) reaction to produce //-mesic nuclei. This 
O I reaction has been used to observe deeply bound pionic states and proven to be 

powerful recently. We calculate r/-mesic bound states in the nucleus using an optical 
OO ■ potential and their formation cross section with the Green function method. Then, 

we carefully check the experimental feasibility. We find that r/-mesic nuclei can be 
observed experimentally using the (d,^He) reaction. We also mention the possibility 
to study the formation of w-mesic nuclei. 
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X 

Bound states of rj mesons in nuclei (j^-mesic nuclei) are interesting objects 
which have not been observed so far experimentally. The //-meson is a mem- 
ber of the SU(3) nonet of pseudoscalar mesons and believed to be one of the 
Nambu-Goldstone bosons of spontaneous chiral symmetry breaking. The ori- 
gin of the mass of the Nambu-Goldstone bosons is studied theoretically in 
terms of a symmetry breaking pattern[l]. Since chiral symmetry is expected 
to be partially restored at finite density[2], it is very interesting to study the 
behavior of mesons, especially their masses, in the nucleus. 

Since the 77-nucleon scattering length is dominated by the s-wave part due to 
the strong coupling to the A^*(1535) resonance, the ry-nucleus optical potential 
has a large s-wave part. Thus by spectroscopic studies of T^-mesic nuclei, we 
can expect to obtain precise information on the s-wave potential, which is 
equivalent to the mass shift of the ?7-meson in the nucleus. We also expect to 
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get new information on the properties of the A^*(1535) resonance in nuclear 
matter by studying the ry-nucleus optical potential. 

The in-medium behavior of tt and K mesons, which are also Nambu-Goldstone 
bosons, has been reasonably well understood from scattering as well as mesic 
atom data. In contrast, there is little experimental information on the in- 
medium properties of the t] meson. Existence of ?7-mcsic nuclei was suggested 
theoretically by Haider and Liu[3]. They systematically investigated ?7-mesic 
nuclear states and proposed to use the (7r+,p) reaction for their formation. 
Experimental attempts to find a bound state in this reaction led to negative 
results[4]. On the other hand, the cross section for 77 meson production in 
d(p,^He)r7 reactions at threshold was found to be largc[5] and was analyzed 
in terms of a quasi-bound 77— •^He system[6]. For the ^He system also the 
existence of a quasi-bound state was suggested [7]. Recent theoretical work 
indicated the existence of ry-mesic nuclei, however their structure is only pre- 
dicted with large uncertainty [8,9]. With the present knowledge, the existence 
of ?7-mesic nuclei is still controversial. 

In this paper we discuss as a new experimental method to produce r^-mesic 
nuclei in (d,^He) transfer reactions on light target nuclei, as recently proposed 
at the GSI heavy-ion synchrotron SIS [10]. The (d,'^He) reaction at recoilless 
kinematics was proven to be a powerful experimental tool by the discovery of 
deeply bound pionic atom formation[ll,12] and to be very useful to extract the 
pion properties at finite density [13, 14]. Here we calculate the structure and for- 
mation cross section of ?7-mesic nuclear states theoretically and investigate the 
experimental feasibility. We should mention here that the same experimental 
technique can be used to produce other mesic nuclei. The u! meson is expected 
to be around 16% lighter at normal nuclear density than in free space[15] and 
expected to form quasi-bound states. The a;-mesic nuclear states can also be 
observed by the (d,^He) reaction in principle as discussed in Ref. [10]. 

In order to study the structure of 77— mesic nuclei, we used the first-order in 
density 77-nucleus optical potential. 



where a^Ar is the 77— nucleon scattering length, jj, is the reduced mass of the 77 
and is ~ ttt.^ for heavy nuclei, M^v is the nucleon mass, and p is the nuclear 
density. 

There exist several recent estimates on the 77A^ scattering length: 




(1) 



a^N = [(0.717 ± 0.030) + i(0.263 ± 0.025)]fm [16], 
= [(0.751 ± 0.043) + i(0.274 ± 0.028)]fm [17], 



(2) 
(3) 
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^ (0.52 + i0.25)fm [7], (4) 
^ (0.20 + i0.26)fm [18]. (5) 

As shown, the first two theoretical estimates agree fairly well with each other. 
The third value was deduced from an experimental study of d(p,^ He)77 and 

d(d,^ He)?] reactions[7]. In all cases, the ?7-nucleus optical potential is expected 
to be attractive. For an illustrative purpose, let us take /i = = 547 MeV, 
Mn = 939MeV and po = 0.17fm-3 and a^jv = 0.717 + 0.263 iim. We then 
obtain 

V{r) = -(86 + 32i)p(r)/po MeV, 
which is indeed strongly attractive. The imaginary part W — —V/2 is appre- 
ciable, but small enough compared with the real part. 

With this potential, we calculated the 77-nucleus binding energies and widths 
for various nuclei in a conventional way of solving the Klein-Gordon equation. 

The vector part of the potential, which in general must be taken into account, 
was ignored in these and the following calculations, based on the assumption 
that the rj — N interaction is dominated by the s— wave component. A Woods- 
Saxon form of the nuclear density profile was used, where nuclear radii and 
diffuseness were taken to be = 1.18^4^/^ — 0.48 fm and a — 0.5 fm, respec- 
tively. The results are shown in Table 1 for the case of a^^j^f = (0.717 + 0.263i) 
fm and in Table 2 for the case of a^Tv = (0.20 + 0.26i) fm. We find that in 
the former case the half widths are comparable or smaller than the binding 
energies and/or level spacings, so that it is justified to interpret these states 
as quasi-stable 77— mesic nuclear bound states. 

Similar to the case of deeply-bound pionic atom production, it is possible to 
produce ?]— mesic nuclei near the recoilless condition using the (d,^He) reaction 
on nuclear targets. This is illustrated in Fig. 1, which shows the momentum 
transfer q vs. the incident deuteron kinetic energy Tj for a typical light tar- 
get nucleus (^Li in this case). The use of recoilless kinematics is essential to 
suppress the quasi-free continuum rj production and to enhance the 77— mesic 
nuclear production signal. The recoil-free r] condition is satisfied at ~ 3.6 
GeV and can be fulfilled at GSI-SIS where the maximum deuteron kinetic 
energy T^"^ is 4 GeV. 

We now estimate the reaction cross section by using the nuclear response 
function S{E): 

\ / A{d,''He)r,(A-l) \ / p((i,-^ife)r; lri,Jn,J 



(, \ lab 
^^^^ ^ is the elementary cross section in the laboratory frame. 

A comprehensive and consistent approach to calculate the response function 



3 



S(E) for a system with a large imaginary potential was formulated by Mori- 
matsu and Yazaki[19]. This method uses the Green function G{E; f, f) defined 
as 

G{E;ry) =< p-^|</.,(f)—-^— -01,(0 >, (7) 



where 0], is the rj creation operator and \p~^ > is a proton hole state. The 
Hamiltonian Hr^ contains the //—nucleus optical potential. Since we used energy- 
independent local potentials in the present calculation, we can obtain a simple 
expression for the Green function as 



G{E;f,f) = Y: yi:,rn,ir)\,rn,ir')Gi^iE-,r,r') (8) 

If) )^77 

Gi^{E- r, r') = -2i^kui^{k, r^)vl+\k, r>), (9) 

where ui^ and Vi^^ respectively are the radial part of the regular and outgoing 
solutions of equation of motion. Using the Green function, the response can 
be calculated as 

S{E) = --Im J2 I d^rda(fr'da'f{f,a)G{E-r,r')f{P,a'). (10) 



We define /(r, a) as 

/(r, a) = x)mimM)K^) ® "i^^Ar. ^)]7MX^(r), (11) 



where Xi ^-^d Xf respectively denote the projectile and the ejectile distorted 

waves, if) is the proton hole wavefunction and ^ is the spin wavefunction intro- 
duced to count possible spin directions of the proton in the target nucleus. The 
numerical values of S{E) were evaluated by using the eikonal approximation 
as in the case of deeply-bound pionic atoms [20]. 

The elementary cross section for 77 production which appears in Eq.(6) can be 
inferred from the energy dependence of the p(d,'^He)?7 cross section measured 
at SATURNE[21] in the d(p,3He) reaction. At Tp = 1.75 GeV (this proton 
kinetic energy corresponds to the recoilless 77 production in the p(d,^He)77 
reaction), the cm. cross section {da/dfl)crn is 3 nb/sr. This can be translated 
to the d + p laboratory- frame cross section via 

da- _ / piab{^Re) y da 

dQlab \Pcm(^He) J dricm' 
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and the approximate elementary cross section was deduced to be 150 nb/sr. 

In Fig. 2, we show the calculated spectra using the Green function method 
described above. The results are shown for the ^Li target (left panel) and for 
the ^^C target (right panel), for different potential parameters. The top and 
middle figures respectively correspond to the rjN scattering lengths of Eq.(2) 
{y{r-) = -(86 + 32i)p{r)/po MeV) and Eq.(4) {V{r) = -(62 + 30i)p{r)/po 
MeV). The bottom figures are for the potential with no binding, V{r) — 
-30ip(r)/po MeV. 

In solid lines, the expected double-differential forward (0°) cross sections are 
shown. The dashed and dash-dotted lines respectively show the contributions 
from the (^3/2)^^ ® {'^P)v ("51/2)^^ ® substitutional configura- 

tions. These two configurations contribute dominantly to the Q-value spectra, 
although we in fact calculated contributions from other partial waves (up to 
/ = 6) and confirmed that there are no significant contributions from partial 
waves beyond / = 6. 

The vertical hues indicate the rj production thresholds; for the ^Li case, the 
threshold is at Qq = -552 MeV while it is at Qq = -558 MeV for the 
case. The 1] binding energy can be deduced from the reaction Q value as 
(for the sake of simplicity we ignore the nuclear recoil energy, which is small 
near the recoilless condition) : 

Q-Qo^ Br, - {Sp{jp) - Sp{p3/2)), (13) 

where {Sp{jp) — Sp{p2,/2)) is the proton hole energy measured from the ground 
state of the residual nuclei. Hence, for the r] states coupled to the (^1/2)^^ 
configuration, the (si/2)p^ — (P3/2)p^ energy differences (14 MeV for '^Li and 
18 MeV for ^^C) taken from ref. [22] was added when calculating these spectra. 

Note that the ground state of the 77— nucleus system for these light p— shell 
targets would have the {p?,/2)p^ ®{^s)n configuration, but this component does 
not contribute to the Q-value spectra near the recoilless condition. Instead, 
the dominant contribution comes from the (^3/2)^^ ® (2p)^ configuration, and 
we can determine the 77— nucleus potential from the location of the 2p peak. 
This (p3/2)p ^ ® {'^P)v component is more dominant in the ^^C case because 
there are four P3/2 protons in a ^^C nucleus as compared to only one in a ^Li 
nucleus. 

In Fig. 3, we show the calculated spectrum for a heavier target, ^°Ca. The po- 
tential strength corresponds to the scattering length in eq. (4). Due to the re- 
coilless condition, dominant contributions come from the substitutional states. 
However, as can be seen in the figure, several proton-hole states contribute to 
a broad maximum in the Q value spectrum which makes the interpretation 



5 



difficult. We tlierefore conclude that light nuclei like those in thep-shell region 
are most suitable to search for bound nuclear states of rj mesons. 

In order to generate realistic Q- value spectra we also have to include back- 
ground contributions from other reactions, which we will discuss in the fol- 
lowing. 

Wc first note that the background due to '^Hc formation without meson pro- 
duction (such as due to coalescence) must be negligible. Composite particle 
production with a few GeV protons incident on nuclear targets was studied by 
Tokushuku et al. [23] at KEK. By extrapolating their results on the deuteron 
spectra in the Al(p,d) reaction at Tp = 3 GeV to ~ 4 GeV/c, we found the 
cross section to be around 10~^[nb/(MeV/c)sr], which is negligibly small. 

The continuum background for 77 production was estimated by using the data 
of Berthet et al. [21]. By relating the number of events for r] production as 

shown in figure 1(b) of Ref. [21] to the tabulated cm. cross section, we esti- 
mated the cm. continuum background level to be ~ 0.09nb/sr/MeV at Tp = 2 
GeV (Trf = 4 GeV). This corresponds to dPa/dEdiliab ~ 4.5 nb/sr/MeV in 
the (d,^He) laboratory frame. 

The d(p,^He)7r"'"7r~ data near the t] threshold by Mayer et a/. [24] show that 
the continuum background due to tt+tt^ production is nearly flat across the rj 
production threshold. We therefore ignored the possible Q-value dependence 
of the continuum background. 

In order to evaluate the continuum background for the case of nuclear tar- 
gets, wc need to calculate the distortion effects of the deuteron and the '^He 
in the target nucleus. For this purpose, we summed up the effective numbers 
for all final state configurations of proton-hole and mesonic states. The calcu- 
lated effective proton numbers had negligible energy dependence. We further 
assumed that the contributions of target protons and neutrons to the back- 
ground is identical, and hence multiplied the calculated effective number by 
a factor A/Z. This total effective number is expected to be a good estima- 
tion of the distortion effects to the projectile and the ejectile. And this is also 
expected to be consistent with the estimation of signal cross sections. The 
effective nucleon number contributing to the background was calculated, and 
was used to estimate the constant background level as: 



A^eff = 0.253 X 7/3 



^ 4 5 X Q X 7/3 = 2.7nb/sr/MeV for^Li 

\ / background 



A^eff = 0.373 X 12/6 
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( 



dEdn 



da 



) 



background 



4.5 X 0.373 X 12/6 = 3.4nb/sr/MeV for^^C. 



In Fig. 4, we show expected Q- value spectra for the "Li case assuming 100 hours 
of beam time at GSI, using the FRS as ^He spectrometer. As shown, we expect 
the peaks to be clearly visible above background, and the spectra are sensitive 
enough to differentiate between various 77— nucleus potential parameters. The 
experimental setup will be similar to the one used for the study of deeply- 
bound pionic atoms. For the estimate we used a target thickness of Ig/cm^, 
a deuteron beam intensity of 3 x 10^°/sec at 3.5 GeV incident energy, and 
= 2.5 X 10~^sr for the FRS acceptance; all these parameters are achievable 



Here, we would like to mention that the (d,^He) reaction is also well-suited 
for the production of a;-mesic nuclei. In Fig. 5, we show calculated spectra for 
uj production in the (d,^He) reaction on ^Li at = 3.8 GeV (also possible at 
GSI), in which we compare three different cj- nucleus optical potentials. At this 
incident energy, however, the recoilless condition is not satisfied, and we hence 
find that the contributions from substitutional states are not dominant, and 
that the quasi-free process makes a large contribution in the unbound region. 
Although the identification of boTind states appears to be difficult, the effect 
of an attractive tu-nucleus potential is noticeable in the bound region of the 
Q value spectrum. 

In order to see the effect of the recoil free condition for u; production we show 
in Fig. 6 the calculated Q value spectrum at 10 GeV incident deuteron energy. 
We assumed the elementary u production cross section to be 450 nb/sr. In 
this case, as in the /^-production spectrum at = 3.8 GeV, the dominant 
contributions result from substitutional states. With less configurations con- 
tributing the a;-nucleus optical potential may be deduced from the spectral 
shape more directly. 

One should note that the (d,^He) reaction is the nucleon pickup reaction, 
allowing for recoil free 77 and cu production, with the lightest projectile and 
ejectile nuclei. The (p,d) reaction which may also satisfy the recoil free con- 
dition at appropriate incident energy, however does not allow to separate the 
ejectiles from the beam particles in a magnetic spectrometer due to the same 
magnetic rigidity in the case of vanishing momentum transfer. 

In conclusion, we find that the proposed recoilless (d,^He) reaction is a promis- 
ing tool to study the 77— nucleus system, and we should be able to determine 
the r]— nucleus potential (and the possible rj mass shift in nuclei) from the 
Q-value spectra. In principle, this method can be extended to study the be- 
havior of other mesons such as a; in nuclei, although there is at present no 



at GSI. 
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facility in the world where one can study uj production near the recoilless con- 
dition (T^ ~ 10 GeV is required). The proposed method is complementary to 
studies of vector mesons in nuclear matter by analysing their invariant mass 
spectrum in the dilepton decay channel, such as — e^e~. A proposed exper- 
imental study of ri and lo production in the (d,^He) reaction at low mometum 
transfer [10] was recently approved at GSI. 

The authors would like to thank H. Toki, T.-S. H. Lee, K. Itahashi, H. Gilg, 
F. Klingl, T. Waas, W. Weise, P. Kienle and T. Yamazaki for helpful dis- 
cussions. This work is supported in part by the Grant-in-Aid for Scientific 
Research, Monbusho, Japan. 
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2 4 6 8 10 

Incident deuteron energy T^j (GeV) 

Fig. 1. The momentum transfer q vs. incident deuteron kinetic energy in the 
Li(d,^He)^He reaction. The three curves respectively correspond to rj binding en- 
ergies of 100, 50 and MeV, as indicated 
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Fig. 2. The calculated rj production spectra for the Li(d,^He) reaction (left) 
and for the ^'^C(d,^He) reaction (right) at Tii = 3.5 GeV, for three different 
rj-nucleus optical potential parameters; (top) V = —(86 + 32i)p/po MeV, (mid- 
dle) V = -{62 + 30i)p/po MeV, (bottom) V = -30ip/po MeV. The vertical lines 
indicate the r] production threshold Q-value (Qq = —552 MeV for the Li case and 
-558 MeV for the C case). In each figure, the contribution from the {0p3/2)p^ ^Pr^ ^'s 
shown in a dashed curve, the (0si/2)p ^ (X" contribution is shown in a dash-dotted 
curve, and the solid curve is the sum of ij-partial waves up to I = 6. The contin- 
uum background contributions are estimated to be about 2.7 nb/sr/MeV for the 
^Li target and 3.4 nb/sr/MeV for the ^"^C target (see text). 
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1.2 



T] production on Ca 
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Fig. 3. The calculated t] production spectrum for the ^^Ca(d,^He) reaction at 
Td = 3.5 GeV, for V = —(62 + 30i)p/po MeV. The labelled curves denote con- 
tributions from the following conhgurations: a) [(1^5/2)^^ '^drj], b) [(1^3/2)^"'^ "Xi^jj], 
c) [(2s)-i (g) Sr,], d) [{lp)p^^Pn] and e) {Is)'^ (g) Srj. 
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Fig. 4. Expected Q value spectrum for the '^Li(d,^He) reaction near the rj produc- 
tion threshold for 100 hours of running at GSI (see text for assumptions of the 
experimental conditions). 
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CO production on ^Li 
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Fig. 5. The calculated oj production spectra for the ^Li(d,^He) reaction at 
Td = 3.8 GeV, for three different iv-nucleus optical potential parameters; (top) 

V = -(100 + 70i)p/po MeV, (middle) V = -(100 + 35i)p/po MeV, (bottom) 

V = —35ip/pQ MeV. The vertical lines indicate the oj production threshold of 
Qo = —787 MeV. In each figure, the contribution from the (0^3/2)^^ "S) Pui is shown 
in a dashed curve, the (0si/2)p^ ® contribution is shown in a dash-dotted curve, 
and the solid curve is the sum of partial waves up to I = 6. The continuum back- 
ground contributions are estimated to be about 7.7 nb/sr/MeV. 
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CD production on ^Li 
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Fig. 6. The calculated uj production spectrum for the Li(d, He) reaction at = 10 
GeV. The elementary cross section was assumed to be 450 nb/sr. The contribution 
from the {0p3/2)p^'^Pv shown in a dashed curve and the (0si/2)p ^®s»7 contribution 
is shown in a dash-dotted curve. The solid curve is the sum of to— partial waves up 
to I = 6. 
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Tabic 1 



A-dependcncc of the rj-nucleus binding energies and widths. We use a^N = 0.717 + 



0.263i [fm] as r]-N scattering length. 
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